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The electron spin coherence in n-doped and undoped, self-assembled CdSe/Zn(S,Se) quantum 
dots has been studied by time-resolved pump-probe Kerr rotation. Long-lived spin coherence per- 
sisting up to 13 ns after spin orientation has been found in the n-doped quantum dots, outlasting 
significantly the lifetimes of charge neutral and negatively charged excitons of 350 — 530 ps. The 
electron spin dephasing time as long as 5.6 ns has been measured in a magnetic field of 0.25 T. 
Hyperfine interaction of resident electrons with a nuclear spin fluctuations is suggested as the main 
limiting factor for the dephasing time. The efficiency of this mechanism in II- VI and III-V quantum 
dots is analyzed. 
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I. INTRODUCTION 

The electron spin coherence in semiconductor quan- 
tum dots (QDs) has been extensively studied experimen- 
tally and theoretically in recent years. [ll-Q This activity 
is driven by prospective applications of the electron spin 
in quantum information processing requiring methods of 
initialization, manipulation, storage, and readout of spin 
coherence. A vital problem in this respect concerns the 
coupling of the electron spin to its environment, deteri- 
orating the spin dynamics. The strength of this interac- 
tion determines the transverse spin relaxation time 72, 
limiting coherent state manipulation. In QDs the elec- 
tron wave function is strongly spatially localized by the 
three-dimensional confinement potential. On one hand 
this leads to discrete energy levels and makes the electron 
spin weakly susceptible to the main spin relaxation mech- 
anisms in bulk semiconductors related to the spin-orbit 
interaction. 3 On the other hand it enhances the Fermi 
contact hyperfine interaction, which therefore is consid- 
ered as main relaxation mechanism for QD electron spins. 
@, @ The hyperfine interaction is controlled by: (1) the 
nuclear isotope composition, (2) the hyperfine coupling 
constants, (ii) the number of lattice nuclei, and (iv) the 
nuclear spin magnitudes. Therefore, a proper choice of 
the QD semiconductor material may minimize this in- 
teraction, potentially prolonging the transverse spin re- 
laxation time. In this respect II- VI semiconductors are 
an interesting alternative to III-V materials. For ex- 
ample, in CdSe the abundance of nuclear isotopes with 
nonzero spin is 7.6% for Se and 25% for Cd only, com- 
pared with the 100% for In, Ga and As in the widely 
studied (In,Ga)As dots.Q Therefore, the electron spin 
dephasing by the hyperfine interaction should be consid- 
erably weaker in CdSe dots. 

Coherent electron spin dynamics was studied by time- 
resolved Faraday rotation for colloidal CdSe QDs fa -fl^l 
and Cd(S,Sc) QDs in semiconductor-doped glasses [l3|. 



both of wurtzite-type. Longitudinal spin relaxation times 
of Ti ^ 20 and transverse spin dephasing times 
Tj* ^ 3 ns were reported. [1, Q The measured Faraday 
rotation traces show complicated quantum beat patterns 
due to the QD-anisotropy typical for wurtzitc structures 
(leading to a complex exciton fine structure) and also 
due to the random dot orientation in the ensemble, so 
that the optical selection rules become undefined. These 
factors complicate the theoretical description of coher- 
ent spin dynamics in colloidal QDs, in addition to the 
influence of surface states. Q 

Some of these complications become resolved in 
epitaxially-grown self-assembled QDs. First, the influ- 
ence of surface states is eliminated. Second, modulation 
doping allows one to fabricate charged QDs with a single 
resident electron per dot. [31 The spin dynamics of resi- 
dent electrons is not limited by radiative exciton recom- 
bination as in empty QDs. Resident electron spin coher- 
ence can be generated via an intermediate trion state, as 
suggested theoretically [11] and demonstrated experimen- 
tally in (In,Ga)As/GaAs QDs. [11, Ull Another advan- 
tage of self-assembled QDs is their well-defined crystallo- 
graphic orientation on the substrate, from which optical 
selection rules are established. Also the crystal structure 
can be controlled by the choice of substrate material, 
e.g., grown on a GaAs substrate, CdSe QDs have cubic 
symmetry instead of wurtzite symmetry, simplifying the 
exciton fine structure. 

In this paper we report on timc-rcsolvcd pump-probe 
Kerr rotation (KR) studies of the electron spin coherence 
in n-doped and undoped self-assembled CdSe/Zn(S,Se) 
QDs with a cubic crystal structure. We observe long- 
lived spin beats related to the coherent precession of 
electron spins about a transverse magnetic field. The 
coherence persists over times exceeding the interval be- 
tween subsequent laser pulses of ~13 ns. Analysis of the 
KR traces as function of magnetic field gives informa- 
tion about dephasing times and g factors of the resident 
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QD electrons. Also the temperature dependence of the 
dephasing time was investigated, and the results are com- 
pared to literature data, especially for colloidal wurtzitc 
CdSe QDs. 



II. EXPERIMENT 

The investigated CdSe/Zn(S,Sc) QD multilayer struc- 
tures were pseudomorphically grown by molecular-beam 
epitaxy on a GaAs:Si buffer layer on top of an n-type 
GaAs (001) substrate. First, a 1.5 //m-thick ra-type Zn- 
MgSSerCl wide-gap confining layer {Eg = 2.96 eV at 
T = 77 K) was deposited to separate the QD region 
from the GaAs buffer. Then, ten periods of CdSe QD 
layers (each having a nominal thickness of 2.5 ML with a 
self-assembled dot density of 5 — 8 x 10^° cm~^) were 
grown, each layer sandwiched between Zn(S,Se)/ZnSc 
short-period superlattice (SL) barriers with a thickness of 
50 nm. The compressive strain induced by the CdSe QDs 
is compensated by the tensile strain intentionally gener- 
ated in the SL barriers, which provides similar growth 
conditions for each QD layer. Finally, the entire struc- 
ture was capped by a 30 nm-thick ZnMgSSe layer. Two 
structures were fabricated: a reference sample with nom- 
inally undoped QDs and a sample with n-doped QDs. 
In the n-doped structure, the middle part of each SL 
barrier with a thickness of 17 nm was doped by CI with 
a concentration resulting in an electron sheet density of 
(3 - 5) X 10" cm-2. We expect that only 10 - 20 % 
of the electrons are transferred to the QDs through the 
undoped 17 nm-thick SL spacers. From these numbers 
we expect that the QD ensemble represents mostly a mix 
of charge-neutral and singly-charged QDs, but there may 
be also a small fraction of doubly or higher charged dots. 

The pulsed laser source used in the experiments was an 
optical parametric oscillator pumped by a mode-locked 
TirSapphire laser. The laser system generates pulses with 
a duration of 1.6 ps at a repetition rate of 75.6 MHz, 
corresponding to 13.2 ns pulse separation. The photon 
energy could be tuned in the range 2.38 — 2.47 eV. The 
samples were held in a helium bath cryostat equipped 
with a superconducting split-coil magnet. Magnetic fields 
B up to 7 T were applied normal to the structure growth 
z-axis (Voigt geometry, B_Lz). The sample temperature 
was varied from 1.6 up to 50 K. 

In the KR experiment the laser beam was split into a 
circularly polarized pump and a linearly polarized probe. 
Both beams were focused on the sample with a spot di- 
ameter of ^150 fim. To avoid dynamic nuclear polariza- 
tion, the helicity of the pump beam was modulated at 
50 kHz rate using a photoelastic modulator. The pump 
pulse excites carriers with spins polarized along the z- 
axis. The subsequent coherent spin dynamics in form 
of Larmor precession about B is measured by the rota- 
tion angle Qkr of the polarization plane of the probe. 
To detect &kr, a- homodyne technique based on phase- 
sensitive balanced detection was used.llJI 



The spectrally-resolved recombination dynamics of the 
optically excited electron-hole pairs could be detected in 
time-resolved photoluminescence (TRPL). For that pur- 
pose the samples were mounted in a helium flow cryostat 
at T = 8 K. The PL emission was dispersed by a 0.5 m 
spectrometer and detected by a synchroscan streak cam- 
era with an S20 photocathode. The overall time resolu- 
tion of this setup was 15 ps. 



III. RESULTS AND DISCUSSION 

The low temperature photoluminescence spectrum of 
the n-doped CdSe/Zn(S,Se) QDs under quasi-resonant 
excitation (see arrow indicating the laser photon energy) 
is shown in Fig. [IJa). The QD emission due to radia- 
tive recombination of negatively charged excitons (trions) 
in singly- charged dots and of charge neutral excitons in 
empty QDs is centered at 2.42 eV. The emission is charac- 
terized by considerable spectral broadening of ^50 meV 
due to QD inhomogeneities. A similar PL spectrum is 
obtained for the undoped QDs (not shown). 

The exciton recombination dynamics was measured by 
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FIG. 1: (a) Photoluminescence spectrum (line) and photo- 
luminescence decay time dispersion (circles) of the n-doped 
CdSe/Zn(S,Se) QDs measured for quasi-resonant excitation 
at 2.46 eV photon energy, (b) Example of the PL intensity 
decay at 2.38 eV after pulsed excitation, plotted on a loga- 
rithmic scale: circles - experiment, line - fit. 
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TRPL. An example of the PL decay at 2.38 eV is shown 
in Fig. [TJb). The decay can be fitted well by a single 
exponential as shown by the solid line with a charac- 
teristic time of tr = 530 ps. The PL decay time varies 
across the emission band: tr increases from 350 ps on the 
high energy flank at 2.43 eV to 530 ps on the low energy 
side at 2.38 eV, see circles in Fig. [TJa). The observed 
systematic increase of the decay times with decreasing 
transition energy is in agreement with previous measure- 
ments on epitaxially-grown CdSe/ZnSe QDs, both for 
ensembles [TtI and single dotsfl^- The increase in the 
exciton radiative lifetime with increasing localization en- 
ergy can be explained by a decrease of the exciton co- 
herence volume in real space and, hence, its increased 
spread in K spacefl^. We note here, that in the stud- 
ied structures the effective barrier energy provided by 
the super lattice minibands is 2.84 eV, which implies the 
strong confinement regime for the confined electron and 
hole wavcf unctions. 

The spin dynamics was investigated by pump-probe 
Kerr rotation. Figure [2]^a) shows the KR trace for the 
71-doped QDs measured at B = 0.5 T. The laser photon 
energy was tuned to 2.44 eV for resonant excitation of the 
trion and exciton states in the QDs. The KR signal shows 
features typical for structures containing singly-charged 
QDs or quantum wells with a low density electron gas.[lj| 
At positive time delays relative to the pump pulse at time 
zero, a periodically oscillating signal with exponentially 
decreasing amplitude for increasing delays is observed. 
Remarkably the signal amplitude is still strong at delays 
of 4.5 ns, i.e., at times strongly exceeding the charac- 
teristic trion lifetimes, that do not exceed 0.6 ns, see 
Fig.[TJa). This allows us to conclude that the KR signal 
is dominated by the contribution of the resident electrons 
in singly-charged QDs, whose lifetime is not limited by 
recombination. A closer look at negative time delays, see 
the red line in Fig. [l{a), reveals spin beats with rather 
small amplitude. 

The spin precession frequencies can be evaluated 
by Fast Fourier Transformation (FFT) as shown in 
Figs, ^h) and M,c)- The precession frequencies for 
negative and positive delays are close to each other: 
UL =8±1 GHz and 7.7±0.2 GHz, respectively. This 
small difference in lol is within the experimental error 
due to the noisy low amplitude signal at negative de- 
lays, which results in a significant broadening of the FFT 
spectrum and a strong noise peak around zero frequency. 
These results also support the incomplete dephasing of 
the electron spin coherence within a laser pulse interval 
of 13.2 ns. 

We consider now the Kerr rotation signals measured on 
the 71-doped CdSe/Zn(S,Se) QDs for different magnetic 
fields in the delay range up to 1 ns, see Fig.[31Ja). At zero 
magnetic field the light induced spin polarization decays 
with three components. Two of them have characteristic 
times of 30 and 300 ps. The third, longest lasting compo- 
nent exceeds the lifetimes of excitons and trions and can 
be assigned to the spin relaxation of the resident elec- 
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FIG. 2: (color online) (a) Time-resolved Kerr rotation trace 
of the n-doped CdSe/Zn(S,Se) QDs vs time delay between 
pump and probe pulses. FFT amplitudes of KR signals an- 
alyzed for negative (b) and positive (c) time delays. Pump 
power 60 W/cm^, probe power 5 W/cm^, laser photon energy 
2.44 eV. 



trons. while the shorter lived components correspond to 
relaxation processes within the trion complexes, which 
arc limited by either spin relaxation or recombination. 
The signals for i? > T are dominated by long-lived 
oscillations with a single frequency corresponding to the 
resident electron Larmor frequency. Also a second non- 
precessing component is present at relatively short delay 
times up to 250 ps. This time range coincides well with 
the trion and exciton lifetimes, so that the non-precessing 
component can be assigned to the spin dynamics of these 
exciton complexes. [20| 

The KR signals from nominally undoped QDs are 
shown in Fig.[3l^b). Their amplitudes are about 20 times 
weaker, but otherwise their appearance is qualitatively 
similar to that of the n-dopcd sample. The signals are 
dominated by long-lived electron spin beats, which can 
be explained by unintentional doping of a fraction of dots 
in the ensemble. No strong hole contribution is seen in 
this sample at short delays. 

In the following we focus on the KR signals of the n- 
doped QD sample and analyze them for delays exceeding 
200 — 300 ps to study the spin coherence of resident elec- 
trons only. For this analysis we used an exponentially 
damped harmonic [31 to fit the KR signals and thereby 
obtain the spin beat frequency and the spin dephasing 
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FIG. 3: KR signals in n-doped (a), and undoped (b) 
CdSe/Zn(S,Se) QDs at difTerent magnetic fields. Pump power 
60 W/cm^ probe power 5 W/cm^, laser photon energy 
2.44 eV. The inset shows the electron g factor dispersion 
across the PL band of the n-doped QD ensemble (circles) 
together with the PL spectrum (line) excited at 2.46 eV. 

time: 

QKRit) = OKRiO) CXp (^-^ ^ COS {uJLt). (1) 

Here Qkb{0) is the KR angle at the moment of pump 
pulse arrival i = ps, T2* is the spin dephasing time, and 
is the Larmor spin precession frequency |2l| . 

First, we discuss the experimental data on the in- 
plain Lande factor of the electron. From the lul 
obtained through the fitting, the Qe can be derived by: 
ge — fujjL/ {^bB). The inset in Fig. |31[a) shows the de- 
pendence of ge on the QD emission energy, g,, has a value 
of about 1.1 with a weak dispersion across the PL band 
originating from the g factor dependence on the band gap 
width [251]. Our experimental technique is not sensitive to 
the sign of Oe , but it is known that g^ is positive in CdSc, 
see Ref. [23] and references therein. 

In Fig. m we compare the ge value for the studied QDs 
with those reported in literature. We chose a presenta- 
tion form in which the ge is given in dependence on the 
energy of optical transition (corresponding to the energy 
gap), as this is the most straightforward way to com- 
pare the experimental results with theoretical predictions 
based on the Roth approach for bulk [2^ . It was shown 
that this approach can be successfully extended also to 
quantum well and quantum dot hetcrostructurcs, where 
the quantum confinement energies have to be added to 
the bulk band gap energy to obtain the optical transition 
energies, see e.g. Refs. P, [23,[25|. In Fig. S] the exper- 
imental data for cubic self-organized QDs are shown by 
the open symbols, and for wurtzite QDs they are given by 
the closed symbols. Our result for CdSc/Zn(S,Se) dots 
(open circle) is in good agreement with the one for cubic 
CdSc/ZnSc dots (open star). [2^ They are also in rea- 
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FIG. 4: Electron g factor as function of optical transition en- 
ergy: n-doped, self-assembled CdSe/Zn(S,Se) QDs (open cir- 
cle), self-assembled CdSe/ZnSe QDs( open star) [2^, wurtzite 
colloidal CdSe QDs (closed circles) and wurtzite Cd(S,Se) 
nanocrystals in glass (closed triangle) [l^. Values for cubic 
bulk materials are shown by arrows: Qe = 0.42 for CdSe and 
1.15 for ZnSe.[22| The dispersion for cubic (Zn,Cd)Se calcu- 
lated after Ref. [22 | is shown by the solid line. 

sonable agreement with the calculations in Ref. [l^j for 
cubic bulk (Zn,Cd)Se, shown by the solid line. These cal- 
culations, in turn, agree well with experimental data for 
(Zn,Cd)Se epilayers grown on GaAs substrates. [22j How- 
ever, the results for the wurtzite CdSe dots shown by the 
closed symbols deviate considerably from that data for 
cubic dots. This difference is due to differences in the 
band parameters for wurtzite and cubic materials, which 
can be also well accounted for in the Roth approach. Q 

The spin dephasing times for the resident electrons 
evaluated from the experimental data in Fig. [3] are pre- 
sented in Fig. [5] as function of magnetic field strength. 
The results for n-doped and undoped QDs coincide well 
and show a strong decrease of the Tg* time from 5.6 ns at 
B = 0.25 T down to 300 ps at B = 7 T. The dephasing 
time is contributed by the coherence time of individual 
spins, T2, and the inhomogeneous dephasing time caused 
by the inhomogencity of the spin ensemble T™'' [l3 |: 



The T™'' time is contributed, e.g., by the nuclear spin 
fluctuations (dominant at very low fields) and by the 
spread of the electron Larmor frequencies (dominant at 
high fields). The latter provides the magnetic field de- 
pendent contribution to the Tj|"'' time, and therefore to 
the time: 

Tj"''(B) r»"'XO) h ' ^ ' 

Here A^e is the spread of the electron g factor for the 
optically excited spin ensemble. One can see from Fig. [5] 
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that the experimental data can be well described by a 
1/B dependence for B > 0.25 T, as shown by the solid 
line assuming a Ag^ ~ 0.0055. This demonstrates that 
the dephasing time for not too small magnetic fields 
is governed by the spread of the electron spin Larmor 
precession frequency. A similar behavior has been re- 
ported for colloidal CdSe QDs.[l,Q Note that this spread 
Age = 0.0055 is less that 0.6 % of the mean value. 

The dephasing time T2 = 5.6 ns at 0.25 T is one of 
the longest reported so far for QDs. It is larger than the 
dephasing times reported for colloidal CdSe QDs of 3 ns 
at B = 0.25 T 0, and than - 2 ns at B = 0.2 T in 
(In,Ga)As/GaAs QDs. [13 

Around zero applied magnetic field the spin dephasing 
in quantum dots is expected to be determined by the 
electron interaction with the nuclear spin fluctuations at 
low temperatures. According to Ref. [5| the electron spin 
dephasing time due to this interaction can be calculated 
according (see Appendix for details): 

T* = V2rA = /i, / ^— . (4) 

Here the sum runs over all types of nuclear isotopes in 
the dots, Nl is the number of nuclei in the QD volume, 
and Ij is the nuclear spin. Aj is the hypcrfine constant, 
which is taken for a unit cell with two nuclei, and yj is the 
probability to find the particular type of nuclear isotope 
(see Appendix for details). Using this equation we can 
evaluate the electron spin dephasing time in CdSe dots. 
To estimate the number of nuclei we used the geometri- 
cal volume of the dot, assumed to be of cylindrical shape 
(QD diameter of 4 — 6 nm and height of 1.4 — 2.1 nm), 
and used hyperfine coupling constants of -37.4 /^eV, - 
39.1 /xcV, and 33.6 /icV for the stable isotopes of ^"'^^Cd, 
^^■^Cd, and ^^Se, respectively, each having spin 1/2 (sec 
Table I and Appendix for details). This gives a Nl in 
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FIG. 5: Spin dephasing time T2 vs transverse magnetic field 
for the n-doped (closed circles) and the undoped (open circles) 
self-assembled CdSe/Zn(S,Se) QDs. 



the range from 630 to 2090 nuclei and leads to a de- 
phasing time of 1 — 2 ns depending on the QD size. For 
better agreement with the measured 5.6 ns we have to as- 
sume that either the QD volume is an order of magnitude 
larger, which can be achieved by assuming, for example, 
a height of 2.1 nm and a diameter of 17 nm or a height 
of 4 nm and a diameter of 12 nm (see Appendix), or the 
interaction with the nuclei is weaker, e.g., due to devi- 
ation of the real situation from the modeled one in the 
exchange-box approximation. 

If the QDs would have the same size, the nuclei- 
induced spin dephasing would be much more efficient in 
GaAs and (In,Ga)As dots compared to CdSe dots, be- 
cause the nuclear spins in III-V dots are considerably 
larger: 9/2 for In and 3/2 for As and Ga compared to 
1/2 for Cd and Se. A comparison of the relevant param- 
eters for CdSe and (In,Ga)As dots is given in Table I. In 
addition, scales inversely with the number of nuclei 
with non-zero spin. The abundance of such isotopes in 
(In,Ga)As is 100% for all nuclear species, while it is only 
25% for Cd and 7.6% for Se. This shortens the dephasing 
time in III-V QDs further. As a result, for dots of the 
same size one would expect the electron spin dephasing 
times to be about 10 times shorter in (In,Ga)As than in 
CdSe dots, see Fig. 7(a). In practice, (In,Ga)As/GaAs 
QDs grown by molecular-beam epitaxy have, however, 
considerably larger sizes than CdSe self-organized dots. 
This reduces the difference in the number of non-zero 
spin nuclei, and therefore also the difference in dephas- 
ing times. 

It is important to mention though that numerical es- 
timates based on Eq. Q give smaller values of T2* than 
experimentally observed for both CdSe and (In,Ga)As 
QDs. Most probably this is related to the limitations of 
the used approach, which is based on the "box model" as- 
suming a constant density of the electron envelope wave 
function at all nuclear sites within the QD. Also possi- 
ble effects of nuclear polarization, e.g. nuclear frequency 
focusing [23, are beyond the model consideration. Ac- 
counting for these effects would require detailed infor- 
mation on the QD shape and substantial computational 
efforts, both of which go beyond the scope of the present 
paper. 

We examined also the temperature effect on the KR 
signal. Figure [5] shows how the dephasing time is chang- 
ing with increasing temperature in a weak magnetic field 
oi B — 0.25 T. The T2 time of the resident electrons de- 
creases monotonically from 5.6 ns at T = 1.6 K (/csT = 
0.14 meV) down to 0.8 ns at T = 44 K (fc^T = 3.8 mcV). 

While a strong temperature dependence of transverse 
electron spin relaxation times has been reported already 
for self-assembled QDs, e.g., (In,Ga)As/GaAs [2^ and 
InP/(In,Ga)P QDs [30|, the underlying mechanisms are 
not trivial and are still not fully understood. They are 
not related to thermal escape of resident electrons, as 
the dot confinement potential considerably exceeds the 
characteristic thermal energies. Also A^e is not expected 
to vary strongly in the scanned temperature range, so 
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FIG. 6: Temperature dependence of the dephasing time T2 for 
the n-doped (closed circles) and the undoped (open circles) 
QDs. Pump power 60 W/cm'^, probe power 5 W/cm^, photon 
energy 2.43 eV. 

that shortening due to related inhomogcneities can be 
most likely excluded. 

Consequently, the shortening of the dephasing time 
T2 with increasing temperature can be related to the 
decrease of T2. Spin-orbit related dephasing is ex- 
pected to be too weak to account for this shortening 
(see Ref. [1^ and references therein). An effective spin 
relaxation mechanism providing electron spin decoher- 
ence in a quantum dot has been proposed in Ref. [sij . 
It suggests phonon-mediated fluctuations in the electron 
spin precession, which are caused by the modulation 



of the longitudinal electron g factor and the hyperfine 
field. Indications to that end were observed also for self- 
assembled (In,Ga)As/GaAs QDs [1^, but more detailed 
studies are required to clarify the role of this mechanism 
for CdSe/Zn(S,Se) QDs. 



IV. CONCLUSIONS 

The coherent spin dynamics of resident electrons in 
self-assembled CdSe/Zn(S,Se) QDs have been investi- 
gated by time-resolved photoluminescence and Kerr ro- 
tation. Long-lived electron spin coherence detectable up 
to 13 ns time delays has been found. The spin dephas- 
ing time of 5.6 ns at i? = 0.25 T is one of the longest 
reported so far for CdSc and (In.Ga)As QDs. The mag- 
netic field dependence of the dephasing time follows well 
a 1/B dependence reflecting a relatively small spread of 
the electron g factor A^e = 0.0055 around the mean value 
of = 1-10 within the QD ensemble excited by the laser 
pulse with a spectral width of about a meV. 
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that for the hypcrfine constants Aj of a particular nuclear 
species different values can be found in literature 0, [32I - 
[34| . The reason for that are differences in defining the 
crystal unit cell for which the Aj are calculated. Further 
the isotope abundance was either included [sj or not [s^l 

The Hamiltonian of hypcfinc Fermi contact interaction 
between electron and nuclear spins is:[35| 



hf 



k 



akiSIk), 



(Al) 



where the sum goes over all nuclei, S is the spin operator 
of the electron, Ik is the spin operator of the k-th. nucleus, 
o-k = 'L>oAi;\ipiRk)\'^ , vq is the unit cell volume, Ak is the 
hypcrfine constant: 



Ak 



16'!T^BfJ'k\UciRk)\'' 

ilk 



(A2) 



with ^.k and Ik being the magnetic moment and spin 
of the fc-th nucleus, the Bohr magneton ^b, the elec- 
tron envelope wave function ')p(Rk) at the fc-th nucleus, 
and Uc{Rk) being the electron Bloch function at the fc- 
th nucleus. We write the electron wave function in the 
form ^{Rk) vo'ip{Rk)uc{Rk), as in Ref. 0. The 
normalization conditions are Jy \Tp{Rk)\'^dv = 1 and 
/^^ |uc(-Rfe)pdf 1. With this definition |uc(-Rfc)P 
I/vq, and therefore one can see from Eq. (A2) that the 
Ak are also ~ 1/fo- 

To estimate the electron spin dephasing time we need 
to calculate the dispersion of the hypcrfine field distribu- 
tion A^. Wc assume that the fluctuating nuclear field 
is sum of the independent contributions of the different 
types of nuclei in the dots. Therefore, the dispersion 
has to be calculated as sum of the dispersions from dif- 
ferent contributions: 



A2 



(A3) 



where the sum runs over the nuclei types. Each individ- 
ual contribution is: 



2v^oI,iI, + l) 



E(A,,^)2|V'(i^,^^)|^ (A4) 



with the sum with ^ running over all nuclei of the same 
type j in a QD. We introduce the dimensionless electron 
density ijj = \uc{Rj)\'^VQ which has its maximum at the 
j-th nucleus, and also the number of nuclei in the volume 
of electron localization Vl ■ 



Appendix A 

Here we give details of the evaluation of the nuclear 
spin fluctuation contribution to the dephasing time for 
an electron spin ensemble in singly charged QDs. Note, 



Nl 



(A5) 



where Vj, = (^J dv\ip(r)\'^^ ^ and n is the number of 
nuclei in the unit cell. For the unit cell with two nu- 
clei (n = 2): Nl = 2Vl/vo, Aj = IGn^iBlijilj /i^Ij^o), 
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where vq is the unit ceh volume. For an unstrained cell 
■Co = ao/4, where Og is the lattice constant. 

One can replace the sum over unit cells by an integra- 
tion and evaluate the sum for all nuclei of the same type 
in a cell: 



1 



3(5eMB)2 Nl 



(A6) 



where yj is the probability to find a nucleus of type j. 
Tjj = XjXj^ where >Cj is the isotope abundance {k S [0, 1] 
as in Ref. (32l |) and Xj is the concentration of the substi- 
tution ions in a ternary alloy, e.g., in Ina;Gai_a;As. For 
binary materials Xj = 1. 

Thereby we obtain an electron spin dephasing time of: 



(A7) 



Here the sum is over the nuclear species. The experi- 
mental value of spin dephasing time corresponds approx- 
imately to the half width at half maximum of the nor- 
mal distribution [2l|. Therefore the relation to the mea- 
sured Tj* can be made through Tj* = v^Ta, where Ta is 
equivalent to the time introduced by Eq.(ll) in Ref. Q. 
The difference in coefficients is accounted for by an other 
choice of the unit cell with 8 atoms in Ref. [5| , instead of 
2 atoms used in our paper. 

The values of rjj for nuclei Ga, As, In, Cd, and Se 
are given in Table I. The data for Ga, As, In, and Cd 
were taken from Refs. [H, [H, Ull. In these papers 
the unit cell containing two nuclei has been chosen, i.e., 
Wo = ao/4. To get 77^ we used ao(GaAs)= 0.565 nm 37| . 
ao(InSb)= 0.648 nm (H, and ao(CdTe)= 0.648 nm |3|. 
For Sc we did not find a value of rysc in literature. We 
assume that the bonds in CdSe arc covalent and take 
?7So = "Hcd for our estimation, as was done in Ref. [33| . 
Values of Aj calculated with unit cells with two atoms 
(with ao(GaAs)= 0.565 nm for In, Ga, As nuclei and 
with ao(CdSc)= 0.6077 nm for Cd, Sc nuclei) are given 
in Table I. 

Numerical calculations give « A:2>.1\/Nl [ps] for 
CdSe (with ao(CdSe)= 0.6077 nm [H). The equiva- 



lent dependence for GaAs is Tj* 
In,Gai_,As: « 6.63^7?!/ ^(T 



6.63-v/A^L [ps] and for 
L2 ~ u.uov-i./ vv--r ll-61a;) [ps]. To 
calculate the last expression we used the oq of GaAs 
(ao = 0.565 nm) also for the ternary alloys, independent 
of xwhich is justified because of the weak dependence of 
ao on In concentration: ao ~ 0.565(1 4- 0.072a;) [nm] [37j . 

Figure [TKa) shows the calculated dependencies of the 
electron spin dephasing time T2 on the number of nu- 
clei Nl in the quantum dot, for CdSe, GaAs and 
Ino.2Gao.8As. Figm-e [Tfb) gives the dependence of T2 
on In concentration in Inj;Gai_2;As for Nl — 10^. 

The estimation of the number of nuclei Nl in CdSe 
QDs (embedded in a ZnSe matrix grown on a (001)- 
orientcd GaAs substrate) is done as follows. We assume 
cylinder-shaped pancake QDs so that volume is given by 
QD Vqd = TTcPh/A, where d and h are diameter and 




. (b) 






ln^Ga^_,As : 
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FIG. 7: (Color online) Model calculations of the nuclei fluc- 
tuation effect on the electron spin dephasing time T^: (a) 
Dependence of T2 on the number of nuclei Nl in the vol- 
ume of electron localization for CdSe, GaAs and Ino.2Gao.8 As; 
(b) Dependence of T2 on In concentration in In^Gai-i^As for 
Nl = 10^ 



height of the QD, respectively. We assume that the 
QD has almost zincblende crystal structure. Each cu- 
bic cell contains 4 Cd-Se molecules, i.e. 8 atoms. Then 



the volume of the cell is vq 
For d ^ A — 6 nm and h 
Nl « 630 - 2090. 



Nl 



d/vo- 



1.4 — 2.1 nm one obtains 



J 
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TABLE I: Parameters of nuclear isotopes used for calculating 
electron spin dephasing times. Only isotopes with nonzero 
spin are given here. All data for nuclear spins, Ij, and mag- 
netic momenta, fj,j , are taken from Ref. [3g| 
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A 1^1 1 Tl H fl T1 PP 

J. \. kj l_l.li.l_J.CLli.\_-C 


77 XlO'^ 


A iipV 


^^^Cd 


1/2 


-0.5943 


0.128 


3.6 " 


-37.4 


^^^Cd 


1/2 


-0.6217 


0.123 


3.6 " 


-39.1 


'''^Se 


1/2 


+0.534 


0.0758 


3.6 


33.6 


"^Ga 


3/2 


+2.016 


0.604 


2.61 " 


38.2 


'^Ga 


3/2 


+2.562 


0.396 


2.61 " 


48.5 


''•^As 


3/2 


+1.439 
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4.42 " 


46 


^^^In 


9/2 


+5.523 


0.0428 


6.35 " 


84.6 


^^''In 


9/2 


+5.534 


0.9572 


6.35 " 


84.8 



"From Ref. p | 
''From Ref. 32 j 
'^From Ref. [sl | 



